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We investigate the Kondo effect in two-dimensional disordered electron systems using a finite- 
temperature quantum Monte Carlo method. Depending on the position of a magnetic impurity, 
the local moment is screened or unscreened by the spin of the conduction electron. The results 
suggest that the Kondo temperature takes different values depending on the position of a magnetic 
impurity. We show that the distribution of the Kondo temperature becomes wide and the weight 
at Tk = becomes large as randomness increases. The average susceptibility shows a weak power- 
law divergence at low temperature in strongly disordered systems, indicating a non- Fermi- liquid 
behavior. We calculate the correlation function between the local moment and the spin of the 
conduction electron. The results are discussed in connection with Kondo screening. 



In low-dimensional disordered systems, the conduction 
electron localizes within a characteristic length called the 
localization length pi El . It was shown that the localiza- 
tion effect enhances the electron correlation effects sig- 
nificantly, leading to conspicuous features particularly in 
low-dimensional systems || f|, ||, ||, 0, ||. The Kondo 
effect is a typical phenomenon caused by the electron 
correlation around a magnetic impurity. The Kondo 
effect in disordered systems was investigated using the 
perturbation expansion in the weakly localized regime 
§, H, 0, |l|. On the basis of the s - d model with 
random potentials, Ohkawa et al. || showed that the 
Kondo logarithmic terms for conductivity and suscepti- 
bility are modified into anomalous terms. Within the 
most divergent approximation, it was shown that higher 
order corrections consist of the product of such anoma- 
lous terms and the Kondo logarithmic terms, and that 
the latter is scaled into the same Kondo temperature as 
that in an ordinary Kondo system pd| ]. The tempera- 
ture dependence of the conductivity and the suscepti- 
bility around and below the Kondo temperature in the 
weakly localized regime was calculated using the 1 IN ex- 
pansion based on the Coqblin-Schrieffer model |Tl] ] and 
the perturbation expansion based on the Anderson model 
p2[ , respectively. 

The Kondo effect in strongly disordered systems was 
studied by Dobrosavljevic et al. by taking account of the 
Coulomb interaction among conduction electrons p~3[ . 
On the basis of a slave-boson approach, they showed 
that, as a result of the localization and interaction ef- 
fects, the Kondo temperature is modified, and that the 
magnetic susceptibility and the coefficient of the linear 
specific heat show divergences as temperature approaches 
zero. In fact, a strong spatial distribution in the Kondo 
temperature was suggested from experimental results of 
strong broadening of the Cu NMR line of UCu5_ x Pd x 
flit . In spite of these findings, the effects of strong ran- 
domness itself on the behavior of a magnetic impurity 
have not yet been fully investigated from a microscopic 
viewpoint. 

In this letter, we study the Kondo effect in two- 
dimensional (2D) strongly disordered electron systems. 



Using a finite-temperature quantum Monte Carlo (QMC) 
method |0|l|, we calculate the susceptibility of a mag- 
netic impurity and the correlation function between the 
local moment and the spin of the conduction electron. 

Let us consider the single-impurity Anderson model 
with on-site random potentials described by the Hamil- 
tonian 



+V 5^(4d>a + H.c.) + Un^n d{ (1) 



where random on-site potentials q are chosen to be a 
flat distribution in the interval [— W, W] under the con- 
dition e i = 0, < i,j > denotes the summation of 
the nearest-neighbor sites, and n da — d\d a . The sys- 
tem is a 41 x 41 square lattice and a magnetic impurity 
is located at i = 0. For W > 3.0, the conduction elec- 
tron is localized with the localization length £(W) < 37.5 
p7[ . Thus, the system for W > 3.0 is probably in the 
strongly localized regime at low temperature. We set 
the condition e d + (1/2)U = 0. In the case without a 
random potential, this condition represents the electron- 
hole symmetry. We use the parameters U = 2.0 and 
V = —1.0 in units of t. The resonance width can be 
estimated as A = irpoV 2 ~ 0.39, where po is the den- 
sity of states at the Fermi energy and can be of the or- 
der of the inverse of the bandwidth (~ 1/8). Therefore, 
the system is in the Kondo regime. The Kondo tem- 
perature without randomness can be estimated as |is[| 
T% = [U 2 C//7r 2 ] 1 /2 e xp[-(^/4)(C//V' 2 - V 2 /U)} ~ 0.14. 

Using the QMC procedure reported in refs. 13 and 14, 
we calculate the susceptibility of a magnetic impurity, 



X(T) 



L-l 



Ar (K*t( T - n dl{n)} K t (t ) - n dL (T )}) , (2) 



where At = (3/L with (3 = l/(k B T) and L is the Trotter 
number. We set gpB/2 = 1. We first numerically diag- 
onalize the Hamiltonian without local correlation effects 



2 



between d electrons to obtain the initial Green functions. 
We then calculate the local correlation effects by the 
QMC procedure, and obtain the Green functions needed 
to calculate physical quantities. We use the Trotter time- 
slice size At = 0.25. An open boundary condition is 
employed. In order to see the effect of the boundary con- 
dition on the low-energy properties, we have carried out 
QMC calculation, applying also the periodic boundary 
condition. The difference between the two results is neg- 
ligibly small. We have run 50000 Monte Carlo sweeps for 
measurements after 1000 sweeps in the warming-up run, 
and have taken an average over the Monte Carlo sweeps. 
Note that the effect of a negative sign problem is negligi- 
bly small because of the conditions ed + (1/2)U = and 

In Fig. 1(a), x{T) versus logT is shown as a function 
of W. A magnetic impurity is located at the center of 
the 41 x 41 square lattice. In the case without random- 
ness, we observe an ordinary Kondo behavior. In this 
realization of a random potential, the susceptibility of a 
magnetic impurity shows a divergence for 1.0 < W < 3.5 
at low temperatures, while it shows a local Fermi-liquid 
behavior for W = 0.5. Note that we have confirmed 
that the data in the latter case can be well fitted by 
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FIG. 1: Susceptibility for < W < 3.5. A magnetic impurity 
is located (a) at the center of the 41 x 41 square lattice, and (b) 
at one of four neighboring sites of the center. Parameters used 
were V = 1 and U = 2 in units of t. The Kondo temperature 
without randomness can be estimated as T% ~ 0.14. 



W_ 

a 



1.0 1.5 2.0 2.5 3.0 3.5 



0.513 0.938 0.981 0.996 1.000 0.999 



TABLE I: Divergence exponents of \ 
1.0 < W < 3.5. 



T~ a in Fig. 1(a) for 



x(T) = aT 2 + b. We evaluate the divergent exponent a 
for x{T) ~ T~ a using data below T£ - 0.14 by the least- 
squares method. The results are summarized in Table I. 
For W = 1.0 and 1.5, the susceptibility of a magnetic im- 
purity diverges with a power weaker than that of a free 
spin, while for W > 2.0, a magnetic impurity behaves 
like a free spin. The results indicate that the impurity 
spin at this site is not completely screened by the spin of 
the conduction electron for W > 1.0. 

Shifting the position of a magnetic impurity twenty- 
four times around the center in the same realization of 
the random potential, we calculate the susceptibility of 
a magnetic impurity. At seven positions, the susceptibil- 
ity shows a local Fermi-liquid behavior irrespective of W , 
while at the remaining seventeen positions, the suscep- 
tibility shows a divergence behavior with increasing W . 
In the former typical example is shown in Fig. 

1(b). In contrast to the results in Fig. 1(a), the suscepti- 
bility at this site shows a local Fermi-liquid behavior for 
< W < 3.5 with monotonic decrease in x(0) as a func- 
tion of IT. In the latter case, a logarithmic divergence can 
be seen at three positions, as well as a weak power-law 
divergence at nine positions. The results for weak power- 
law and logarithmic divergences are depicted in Fig. 2 
for W — 3.5. The divergent exponents are evaluated in 
the same way as in Fig. 1(a) and are presented in Table 



#1 • #2 - #3 " #4 □ #5 
#6 » #7 ' #8 o #9 




FIG. 2: logx(T) versus logT for weak power-law divergence 
for W — 3.5. A magnetic impurity in each result is located at 
the different position around the center. Parameters and the 
random potential realized were the same as those in Fig. 1. 
Solid, broken and dotted lines are fitted by the least-squares 
method. Inset: x(T) versus logT for logarithmic divergence. 
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#1 #2 #3 #4 #5 #6 #7 #8 #9 



0.952 0.938 0.807 0.751 0.597 0.485 0.485 0.414 0.343 



TABLE II: 
2. 



Weak divergence exponents of x ~ "in Fig. 



II. In Table III, the behavior of x an d £i for = 3.5 
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log-D 
0.094 


#6 
-0.454 


log-D 
-0.099 


#9 
0.454 


#2 
2.985 



TABLE III: Behavior of \ an d a for W = 3.5 at each shifted 
position of a magnetic impurity. The average e; is zero. The 
numbers in the top row and the column of the left-hand side 
represent the coordinates along the horizontal and vertical 
axes, respectively. The position (21, 21) is the center of the 
system. FL, log-D and FS denote local Fermi-liquid behavior, 
logarithmic divergence and behavior of a free spin, respec- 
tively. #1 ~ #9 correspond to the results in Fig. 2. 

at each position of a magnetic impurity is summarized. 
From the results, we conclude that the local moment in 
2D disordered electron systems is screened or unscreened 
by the spin of the conduction electron, depending on the 
position of a magnetic impurity. Since y(0j = (2ttTk)~ 1 
with Tk being the Kondo temperature |L8| , this conclu- 
sion indicates that the Kondo temperature has a spatial 
distribution down to Tk = 0. 

On the basis of the results for the susceptibility of a 
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FIG. 3: Distribution of the Kondo temperature for W = 1.0 
and 3.5. The position of a magnetic impurity is shifted 
twenty-four times around the center. T% is the Kondo tem- 
perature without randomness. 




FIG. 4: Average susceptibility for W = 3.5 in 0.05 < T < 
Tk ~ 0.14. The broken line is fitted by the least-squares 
method: Xav = 0.867 T" " 



-0.744 



magnetic impurity, we estimate the distribution of the 
Kondo temperature. When the susceptibility shows a 
Fermi-liquid behavior, we extrapolate x(0) by the least- 
squares method using x(T) = aT 2 + b. The results are 
summarized in Fig. 3. At W = 1.0, the Kondo tem- 
perature is distributed around Tk = 0.9 T K , where the 
largest weight exists, and a small weight lies at Tk = 0. 
At W = 3.5, the distribution of Tk becomes wider and 
the weight at Tk — increases considerably with reduced 
weights for Tk ^ 0. We have shown that the spatial dis- 
tribution of the Kondo temperature can be caused only 
by the effects of a random potential. 

In the disordered system, we take an ensemble average 
over the distribution of random potentials in order to ex- 
tract the universal properties. Judging from our numeri- 
cal results, it seems that the susceptibility of a magnetic 
impurity does not show the universal behavior at low 
temperature. Therefore, in 2D strongly disordered elec- 
tron systems with dilute magnetic impurities, where each 
magnetic impurity acts as a single magnetic impurity, the 
observable susceptibility at low temperature may be ob- 
tained as the arithmetic mean over the susceptibility of 
each magnetic impurity. The average susceptibility thus 
obtained in the random potential for W — 3.5 realized 
in Fig. 1 is shown in Fig. 4. The average susceptibility 
for 0.05 < T < T K shows a weak power-law divergence 
Xav ~ T~ a with a ~ 0.744, indicating a non-Fermi-liquid 
behavior. 

To obtain local information about Kondo screening, we 
calculate the correlation function between the local mo- 
ment at the site and the spin of the conduction electron 
at the i site, 



<S; 



0°i 



>= \< t" d T ~ n d i][n^ - nn] >> ( 3 ) 



where = c| CT Ci CT . In Figs. 5(a) and 5(b), we show the 
spin correlation functions along the vertical and horizon- 
tal axes of the square lattice, respectively. The realized 
random potential and the position of the local moment 
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— — W=0.0, same position as for Fig. 1(a) 
— •— W=3.0, same position as for Fig. 1(a) 
■ - a- - w=3.0, same position as for Fig. 1(b) 
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FIG. 5: Correlation function < SqS* > along the (a) ver- 
tical and (b) horizontal axes of the square lattice. Random 
potentials and the position of the local moment are the same 
as those used in Figs. 1(a) and 1(b). Parameters used were 
V = 1, U = 2, W = 3.0, and = 20. Inset: < S§sf > in a 
wide scope. 



are the same as those used in Figs. 1(a) and 1(b). Pa- 
rameters used were W = 3.0 and (3 = 20. Thus, the 
results shown in Fig. 1(a) (open circles) correspond to 
the unscreened case, while those shown in Fig. 1(b) (full 
triangles) correspond to the screened case. In the former 
case, the antiferromagnetic correlation at the magnetic 
impurity is suppressed due to randomness. The suppres- 
sion of the antiferromagnetic correlation at the magnetic 
impurity is commonly seen in other locations of a mag- 
netic impurity where the susceptibility shows a diver- 
gence. This suppression probably causes the incomplete 
screening. Although the antiferromagnetic correlation at 
the magnetic impurity is not suppressed in the screened 
case, the behavior of the spin correlation function around 
a magnetic impurity is different from that for W = 0; the 
ferromagnetic correlation is reduced. 

In summary, we have investigated the Kondo ef- 
fect in 2D disordered electron systems using a finite- 
temperature QMC method. We have shown that de- 
pending on the position of a magnetic impurity, its lo- 
cal moment can be screened or unscreened by the spin 
of the conduction electron. In the latter case, the anti- 
ferromagnetic spin correlation between the local moment 
and the conduction electron is suppressed around a mag- 
netic impurity. We have demonstrated that the effects of 
randomness induce a distribution of the Kondo tempera- 
ture down to Tk = 0, and that the average susceptibility 
diverges with a weak power law at low temperature in 
strongly disordered systems. 
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